The role of light ellipticity in ionization of atoms by intense few-cycles laser pulses 
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We provide theoretical investigations of the response of the Ar and H atoms to an intense el- 
liptically polarized few-cycle laser pulse, as a function of light ellipticity. The time-dependent 
Schrodinger equation describing the least-bound electron is solved numerically, and differential 
quantities such as the momentum distribution, the electron density in the continuum, and the 
above-threshold ionization spectra are computed. These quantities provide insight into the ioniza- 
tion dynamics and the electron rescattering process as a function of light ellipticity, and reveal great 
similarities between the response of Ar and H to the applied external field. 

PACS numbers: 32.80.Rm 



I. INTRODUCTION 

Experiments are now emerging on strong-field ioniza- 
tion of atoms 1] and molecules by femtosecond, 
elliptically polarized lasers. For an elliptically polarized 
laser, the field has two perpendicular components, and 
the electron dynamics and the resultant momentum dis- 
tributions depend on the relative magnitudes of the field 
components (light ellipticity), pulse duration (few-cycle 
vs. many-cycle ionization regimes), and carrier-envelope 
phase (CEP), see, for example, Refs. [5-7] and references 
therein. 

In a recent experiment [l[ on strong-field ionization of 
the He atom by a few-cycle elliptically polarized laser, 
the momentum distributions in the polarization plane 
of the external field are rotated slightly relative to the 
prediction of the simple-man's model, in which the fi- 
nal momentum of a continuum electron born at time t is 
defined (in atomic units, which are used throughout) as 
kf = - J* t °° E(t')dt' = -A(t) where E(t) is the electric 

field and A(t) the vector potential at time t (cf. Eq. [I}. 
For a few-cycle pulse of duration T, electron emission 
at the peak of the pulse (T/2) dominates the ionization 
process, and the final momentum is thus kf — —A(T/2). 
In Q , calculations of the response of the H atom to a few- 
cycle circularly polarized laser pulse, based on solution of 
the time-dependent Schrodinger equation (TDSE), pro- 
duced a similar angular shift. The angular shift has been 
attributed to the interplay between the Coulombic po- 
tential and the external field Q. 

In the present work, we investigate the response of 
many-electron atoms to intense few-cycle elliptically po- 
larized laser pulses. We consider Ar because its ionization 
potential is comparable to that for the H atom, and this 
would facilitate comparisons with the recent TDSE re- 
sults for the latter [8[ . Note that H and Ar produce very 
similar momentum distributions when probed by a lin- 
early polarized few-cycle laser pulse [|[ , because the ion- 
ization process is not sensitive to the short-range poten- 
tial (the precise electronic structure of the probed atom) . 
Here, we extend the investigation to elliptically polarized 
fields, with light ellipticity [e in Eq. (p}] ranging from 7r/2 



(circular polarization) to 7r/10 (approaching the linear 
polarization limit). Note that whereas circular polariza- 
tion switches off electron rescattering, electron rescatter- 
ing effects are very important in the linear polarization 
limit Here, we discuss the effect of light ellipticity 

on the momentum distributions. Our analysis is assisted 
by above-threshold ionization (ATI) spectra and electron 
density in the continuum. These quantities help us de- 
termine at which e value electron rescattering becomes 
important, and this has implications on the application 
of elliptically polarized lasers to investigations of orbital 
structure and symmetry [H, fll| . 

The following are the main findings of this work. The 
momentum distributions of Ar and H show great simi- 
larities between the responses of the two systems to the 
external field, and to change of light ellipticity. For both 
atoms, the momentum distributions are rotated relative 
to the predictions of the simple-man's model. The an- 
gular shift does not depend very much on the probed 
system, indicating that in the present regime it is deter- 
mined mainly by the long-range part of the atomic poten- 
tial. The angular shift is very sensitive to light ellipticity, 
and it is smallest for the circularly polarized light. By 
investigating the electron density in the continuum and 
the ATI spectra as a function of light ellipticity, we find 
that electron rescattering becomes important at e < 7r/4. 

The paper is organized as follows. The computational 
details are given in Sec. UH the results and discussion in 
Sec. IIII1 and summary and conclusions in Sec. IIV1 



II. COMPUTATIONAL DETAILS 

The external field is defined as E(t) — ~d t A(t), where 
A(t) is the vector potential [defined in the xy— plane, cf. 
Eq. (JTJ)]. For an elliptically polarized light, we write the 
vector potential as 



A(t) = A f(t) 



cos(wi - 
sin(o;i 







lcos(e/2) \ 
lsin(e/2) , 



(1) 



where Aq is the amplitude, u the carrier frequency, <j> the 
CEP value, e the light ellipticity, and f(t) = sin 2 (wt/2N) 
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the envelope for an A^-cycle pulse. In the present study, 
uj = 0.057 a.u. (800 nm wavelength), <fr = —tt/2, and 
N = 3 optical cycles. 

The wave function is expressed in spherical harmonics 

as 

*(»■.*)= E E ^YlrnW, (2) 
1=0 m=-i 

and the TDSE is solved in the velocity gauge with a 
grid representation for the reduced radial wave functions 
fim(r) [HI- The SAE potential describing Ar is taken 
from Ref. [13j . We use an equidistant grid with 4096 
points that extends up to 400 a.u. for Ar (300 a.u. for 
H). Truncating the expansion in <j2j) at l m ax — 40 for Ar 
(35 for H) produces converged results at laser intensity 
of 1.06xl0 14 W/cm 2 . We note that for each I, the az- 
imuthal quantum number runs through 21+1 values mak- 
ing the calculations for an elliptically polarized light fully 
3D and much more time demanding than in the case of 
linearly polarized light. 

The calculations are carried out at e = tt/2, tt/3, 7r/4, 
and 7r/10. The momentum distributions are computed in 
the polarization plane (dP/dk with 9k fixed at tt/2) by 
projecting on scattering states [lij]. Since the external 
field (JT|) is defined in the xy— plane, we only consider con- 
tributions from the Ar 3p2; and 3p y states. For symmetry 
reasons, the 3p z orbital will not contribute to ionization 
in the xy— plane. 

III. RESULTS AND DISCUSSION 

In Fig. [TJ we provide parametric plots of the vector 
potentials (top panels) and momentum distributions in 
the xy— plane (dP/dk at k z = ) for Ar (middle pan- 
els) and H (bottom panels) in elliptically polarized fields 
with fixed intensity of 1.06xl0 14 W/cm 2 and light ellip- 
ticity e = 7r/2, 7r/3, 7r/4, and 7r/10, from left to right, 
respectively. We set dP/dk = for fc 2 +fc 2 < 0.01 a.u., 
in order to achieve a better graphical display. The mo- 
mentum distributions for Ar are incoherent sums of the 
individual contributions from the 3p x and 3p y states. At 
e = tt/2, 7r/3, and 7r/4, both states contribute to the 
momentum distributions. At e = 7r/10, by contrast, the 
major polarization axis of the external field is aligned 
with the x— axis as can be seen in Fig. [TJ and the con- 
tribution from the 3p x state dominates the momentum 
distribution. 

Starting with Ar at e = tt/2, the momentum distri- 
bution shows a ring-like structure with large emission 
probability in a finite momentum region [ky ss 0.75 a.u.], 
characteristics of a short pulse duration |6J. The mo- 
mentum distribution shows a spiral structure that devel- 
ops outward with a counter clockwise rotation, following 
the time evolution of the vector potential (JXJ) - Moreover, 
the momentum distribution is characterized by an angu- 
lar shift, relative to the predictions of the final momen- 



tum based on the simple-man's model (kt — —A(T/2)). 
As mentioned earlier, this shift is due to the interplay 
between the atomic potential and the external field Q. 
From the momentum distribution, one can see that the 
probability of finding electrons with low momenta (0.1 
< kf < 0.5 a.u.) is very small, indicating that the contin- 
uum electron is progressively driven away from the core 
by the external field, and that the electron rescattering 
channel is closed. 

In fact, it has recently been shown by both experiment 
and calculations that investigations of oriented samples 
(atomic [ll[ and molecular Q) by a circularly polarized 
light (for which electron rescattering is negligible) pro- 
vide a unique probe of the orbital symmetry— the orbital 
angular nodes are preserved in the momentum distribu- 
tions. 

For Ar at e = 7r/3, the momentum distribution is still 
characterized by a large emission probability in a finite 
momentum region. However, at a progressively lower mo- 
mentum compared to the calculations at e = tt/2. This 
latter feature is understood based on the simple-man's 
prediction of the final momentum. One can also see that 
the probability of finding low-momentum electron is still 
small, compared to the main emission probability. We 
note that the angular shift of the momentum distribu- 
tion increases by going to smaller e values. Our TDSE 
results show that the angular shifts for Ar are 10° at 
e = tt/2 and 27° at e = tt/3. 

Turning to the e = tt/A case, the regular emission pat- 
tern observed at e = tt/2 and tt/3 is now distorted by 
rescattered electrons and bears a clear signature of low- 
energy electrons typical for tunneling ionization by the 
linear component of the external field. Because of these 
new features, it is no longer meaningful to discuss the 
angular shift at e — 7r/4. At e = 7r/10, the momentum 
distribution is mainly along the k x — axis. However, it is 
not inversion symmetric across k y . Moreover, the radial 
jets and the interference structure typical for ionization 
by a linearly polarized laser pulse [9j are washed out by 
the circular component of the external field. 

The momentum distributions computed for H (bottom 
panels of Fig. [JJ) are generally similar to those for Ar, 
and show the same response to change of light ellipticity. 
For instance, the angular shifts computed for H (9° at 
e = tt/2 and 27° at e = 7r/3) are similar to those for Ar, 
indicating that the shift is not sensitive to the short-range 
potential of the probed atom. There are, nevertheless, 
quantitative differences between Ar and H such as the 
lack of spiral structure in the momentum distribution of 
H at e = tt/2. This is due to the different ionization po- 
tentials (0.5 a.u. for H; 0.57 a.u. for Ar). For a circularly 
polarized light with laser intensity of 1.06 x 10 14 W/cm 2 , 
the Keldysh parameter [l5[ value for Ar (H) is 1.6 (1.5). 
The spiral structure has been reported for H at laser in- 
tensity 5 x 10 13 W/cm 2 corresponding to a Keldysh 
parameter value of 2.1. 

The ATI spectra, obtained by integrating the 3D mo- 
mentum distributions over the angular variables, are 
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FIG. 1: (Color online) 2D momentum distribution for Ar (middle panels) and H (bottom panels) in an elliptically polarized 
field (vector potential in top panels) with e = it/2, 7r/3, 7r/4, and 7r/10 from left to right, respectively. The laser pulse contains 
3 cycles and has a peak intensity 1.06x 10 14 W/cm 2 and frequency 0.057 a.u. In the rightmost panels, the dP/dk scale is 
logarithmic. 



shown for Ar and H in Fig. [21 and provide further support 
for our interpretations of the momentum distributions in 
Fig. [TJ For instance, for Ar and H at e = it/2 (and 
still at 7r/3), the ATI spectra show essentially a single- 
peak distribution, in agreement with the predictions of 
the ionization models [16H18J. At e = 7r/4, by contrast, 
the ATI spectra show signature of low-energy electrons, 
characteristics of tunneling ionization by the linear com- 
ponent of the external held. By comparing with the ATI 
spectrum obtained at e = 7r/2 and focusing on the high- 
energy tail, we find that while electron rescattering effects 
are minimal at e = 7r/3, they enhance the high energy tail 
of the ATI spectrum at e = tt/4. At e = tt/10, the ATI 
peak develops a high-energy shoulder due to rescattered 
electrons. 

The electron density in the continuum is shown for Ar 
and H in Fig. [3l just after the end of the laser pulse. 
At e = it/ 2 and 7r/3, the electron density of the contin- 
uum electron is driven away from the core. By contrast, 
at e — 7r/4, the electron density is smeared out due to 
rescattering. At e — tt/10 the electron density is focused 
along the major polarization axis (x— axis) of the exter- 
nal field. 

From the analysis of the ATI spectra (Fig. [5]) and elec- 
tron density in the continuum (Fig. [3]), we learn that 



electron rescattering is essentially absent for elliptically 
polarized laser pulses with it/ '3 < e < it/ 2. These find- 
ings are important for guiding future experiments on the 
investigation of orbital structure by employing elliptically 
polarized laser pulses, which would benefit from a clean 
orbital signal— without the contribution from rescattered 
electrons. 

IV. SUMMARY AND CONCLUSIONS 

To summarize, we investigated ionization of atomic 
targets (H and Ar) by intense 3-cycle elliptically polar- 
ized laser pulses with light ellipticity (e) ranging from 
the circular polarization case (e = 7r/2) and approaching 
the linear polarization case (e = 7r/10). The momentum 
distributions, ATI spectra, and electron density in the 
continuum are studied as a function of light ellipticity, 
and reveal great similarities between the response of the 
two systems to the external field. Both systems possess 
an angular shift in the momentum distributions, relative 
to the predictions of the simple-man's model, and the 
magnitude of the shift is dependent on light ellipticity. 
From analysis of the ATI spectra and electron density 
in the continuum, we obtain further insight into the ion- 
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FIG. 2: ATI spectra (dP/dE in a.u.) for Ar (top panels) and H (bottom panels) in an elliptically polarized field with e — n/2, 
7r/3, tt/4, and 7r/10 from left to right, respectively. The dashed curves in the top panels denote the ATI spectrum computed 
at e — 7r/2, shifted in energy and renormalized. See caption of Fig. [T]for details regarding laser parameters. 




FIG. 3: (Color online) Electron density in the xy— plane for Ar (top panels) and H (bottom panels) just after the end of an 
elliptically polarized pulse with e = tt/2, 7r/3, 7r/4, and n/10 from left to right, respectively. See caption of Fig. [1] for details 
regarding laser parameters. In the rightmost panels, the density scale is logarithmic. 



ization dynamics and the ellipticity onset of the electron rescattering process. 
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